BOUND STATES IN MONOPOLES: SOURCES FOR UHECR ? 
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Bound states in monopoles are studied through a simplified, Witten-like model. As the overall 
structure is determined in full details, it is shown in particular that only those states having a 
vanishing angular momentum are allowed; for those, the energy spectrum is derived numerically 
and an approximation is set up that allows an easy description in terms of wavefunctions, usefull for 
further applications. The monopoles are then proposed as candidates of ultra high energy cosmic 
rays, an hypothesis that should soon be testable through the Pierre Auger Observatory. 
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I. INTRODUCTION 

Among the various topological defects that may have 
formed during phase transitions in the early universe [Q , 
monopoles are absolutely unavoidable in Grand Unified 
Theories (GUT), according to which there was a time 
when all symmetries were unified in a single semi-simple 
gauge group, since electromagnetism is an exact (unbro- 
ken) symmetry still now. They are however usually con- 
sidered as a sufficient nuisance Q (in the sense that they 
would very rapidely come to dominate the universe) that 
inflation is generally invoked right after their time of for- 
mation so that they are completely diluted H and thus 
unobservable. These results are based on the by-now very 
standard view that they originated at the GUT phase 
transition, i.e., at an energy around 10 15 — 10 16 GeV. 
Many ways out however have been proposed other than 
inflationnary scenarios, among which the Langacker-Pi 
mechanism Q which relies upon using cosmic strings to 
connect the monopoles and anti-monopoles pairs, thereby 
effectively enhancing considerably the decay probability, 
thus reducing the remnant monopole density. Although 
the model, in its original presentation, suffers from many 
drawbacks, it has at least the advantage of proposing an- 
other solution, not involving inflation and letting open 
other alternative possibilities. Domain walls have also 
been used to sweep them away ||; in all cases, unde- 
sirable monopoles are gotten rid of by means of higher 
dimensional topological defects Q . 

Yet another alternative possible solution to the 
monopole excess problem is the simplest one, although 
completely overlooked until recently: it consists in not- 
ing that the monopole density is in fact proportional to 
the fourth power of the energy scale ij at which the sym- 
metry breaking during which monopoles were generated 
took place, given the monopole mass m M is essentially r\ 
times the inverse of the corresponding coupling constant, 



i.e., ~ 137 in the case of electromagnetism (as ought to 
be the case): 
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(here we note Q M the monopole density in units of the 
critical closure density, and h the Hubble constant in 
units of 100 km • s — 1 - Mpc -1 ) so that it suffices to lower 
7] to < 10 9 Gev in order to cure the density problem. 
This is the solution we shall adopt here, for it might also 
provide a useful explanation for the high energy cosmic 
ray mystery Q. (Note in this regard that one can also 
simply assume the monopole overdensity problem to be 
cured somehow and investigate anyway the possibility 
that they give rise to high energy cosmic rays ||].) This 
way, the possibility that monopoles are still present in 
the universe is still reasonnable. 

Once the mass scale is fixed, the relevant physics 
still needs be properly clarified in order to enable one 
to study the interactions between monopoles and other 
particles. In this regard, models have been suggested 
where bound states of scalar |J] , gauge vector bosons |l(]] 
or fermions lO] can form in SO(3) t'Hooft-Polyakov 
monopoles ]12|. Those modify the scattering solutions 
and can enhance greatly the cross sections by means 
of the Callan-Rubakov effect It is the same kind 

of approach we wish to present here, although using 
a simplified Witten like model (generally used to de- 
scribe current-carrying cosmic strings fl^ , |l5|| or domain 
walls jl6|]) and a completely different method for calcu- 
lating the bound state energy levels. Our model, con- 
trary to the other proposed models having bound states, 
presents the advantage (apart from its simplicity) that 
the bound states are present for dynamical reasons, i.e., 
they do not exist for all values of the underlying parame- 
ters. This means that we can consider that particles may 
be trapped in the monopole even though they might be- 
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long to representations that are not directly related to 
that of the monopole fields. Besides, we shall work using 
the Julia-Zee fully numerical solution [[l7| and not the ex- 
plicit BogomoPnyi-Prasad-Sommerfield (BPS) analytic 
limit @. 

In what follows, we shall first define our model and 
show that bound states may be present in monopoles, 
provided their angular momentum vanish, exactly as in 
the Callan-Rubakov mechanism JD^. This not-quite- 
obvious result follows from the requirement that bound 
states do exist in the monopole, which in turn seriously 
constrains the possible couplings between the fields in- 
volved. 

Having explored the full microscopic structure of a 
monopole having bound states, we go on to investigate 
the possibility that, in a fashion similar to Witten su- 
perconducting cosmic string models [p^j , through the ex- 
istence of vortons (rotating superconducting cosmic 
string loops configurations), these monopoles could be 
the source of Ultra High Energy Cosmic Rays (UHECR). 
Their acceleration, once they are coupled to electromag- 
netism, poses no particular problem as long as an exter- 
nal magnetic field exists in the accelerating region un- 
der consideration. Similarly to the model based on vor- 
tons |2(| , they can propagate along huge (cosmologically 
speaking) distances, and thus have the ability to reach 
us fairly easily. Their interaction cross-section with air 
nuclei can be obtained through the Callan-Rubakov ef- 
fet Jl3[ ] and is thus evaluated to be hadronic in nature. 
The expected properties of the resulting UHECR distri- 
bution are essentially those derived in the vorton case ]2C| ] 
and can therefore reproduce the existing data. More 
data, thanks for instance to the Pierre Auger Observa- 
tory project will give a definite answer concerning 
this possibility. 



II. THE MONOPOLE STATE 

The model we shall use in what follows involves the 
symmetry breaking of an SO (3) invariance by means of a 
Higgs field $ belonging to the 3 representation of SO (3), 
coupled to a complex scalar field E which we assume, 
for the sake of simplicity, not to be coupled to the gauge 
field A p of SO (3). This assumption should of course be 
modified when one wants to evaluate the long-range elec- 
tromagnetic interaction of the resulting monopole with 
other particle, but we shall show latter on how this can 
be achieved. With a metric convention having positive 
signature, we have the model 

C = - l -{D^ a )\D^ a ) - i(0 p E)*(0"E) - 

Xrh ,~ ' ' ' 2\2 



($ a <S> a -r) 2 y -/($ a $ a -^)|E| 

-^|E| 2 -^|E| 4 , 
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where the covariant derivative is 

D^a = d^ a - qe ab c A b ^ c , (3) 

and the gauge field strength 

F^^d^-d^ + qe^AlAl. (4) 

A static monopole configuration then has the 
form (ljjl^], in spherical coordinates x l = (r, 6, (f) 



$ a = Mr) — , 
r 
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1 qr 2 13 u 



(5) 
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The field equations for the configuration (||) and (|6|) , with 
the E field not taken into account, i.e. for the ordinary 
't Hooft-Poyakov monopole, read 



1 d , ndh. 2 



ir 1 — ) = ^hK z + X^h(h z - 1) (7) 
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= q 2 r] 2 Kh 2 + -^K(K 2 - 1), 

with boundary conditions 

h(0) = K(oo) = 0, h(oo) = K(0) = 1. 



(8) 



(9) 
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FIG. 1. Amplitude of the field functions h(p) and K(p) vs 

1/2 

the rescalled unit of length p — r\. Curves are shown for 
q 2 /Xrf, — 10~ 2 (solid line), 0.1 (dashed line), 1 (long dashed 
line) and 10 (dot-dashed line). 
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These equations have been solved numerically for vari- 
ous values of the only relevant parameter q 2 /A^ and pro- 
duce the characteristic curves on Figure [|, where h and 
K are shown as functions of the rescalled unit of length 
P = r^/KpV- 

The total energy of the monopole is expressible simply 
in terms of q 2 /\<p and the Higgs field mass m H — ^/X^i] 
as 



E.. = 47rm 
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which is shown on Figure g as a function of q 2 /\$. 
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which will admit bound state eigensolutions (to 2 < m 2 ) 
provided that the potential V satisfies either 



V < 



37? e [0,oo[ ; 



= 



<2/ry 2 



(14) 



dV 
dr 



and V(R) < 0. (15) 



In the latter case, denoting by a prime a derivative with 
respect to the rescalled distance p, one finds that the 
minimum of the potential would be for p such that 



A</, 

V 



1(1 + l) = hh'p 3 



(16) 



On Figure || is plotted the right hand side of this relation 
against p for three orders of magnitude of the relevant 
underlying parameter q 2 /X c f > , and it is seen that in general 
this function is of order unity throughout its range of 
variation. Hence, in order for the condition (16|) to be 
satisfied, it is necessary that the left hand side be also 
of order unity. Assuming A^> and / to have comparable 
values (they are both quartic interaction term coupling 
constants), this leads to the constraint 
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(17) 
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which shows that no bound state is expected for large 
angular momentum. For this reason, we restricted our 
analysis to vanishing angular momentum states £ = 0. 
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FIG. 2. Variation of the total monopole energy as a func- 
tion of g 2 /A^>. 



III. THE SCALAR CONDENSATE 

In this background, we now investigate the behaviour 
of the bosonic field S by first looking at the field equation 
in which the nonlinar term is omitted. Separating space 
and time variables in the form 

E(a^) = CT(r)e fe %„(M), (H) 

the field equation for £ in the monopole background gives 
the Schrodinger-like eigenvalue equation for the ampli- 
tude of this field as 



Act + V(r)a = uj 2 a, 



where 



V(r) 



+ 2/? 7 2 (/ i 2 -l)+m 2 ] 



(12) 



(13) 
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FIG. 3. Value of the function p s h(p)dh/dp which de- 
termines the possibility of bound states for nonzero val- 
ues of the angular momentum i vs the rescalled unit of 
length p = r\^ 2 r/. As on Figure |l], curves are shown for 
q j\<t> — 10~ 2 (solid line), 0.1 (dashed line), 1 (long dashed 
line) and 10 (dot-dashed line). 
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Yet another way to convince oneself that the bound 
states should be restricted to those with vanishing an- 
gular momentum consists in investigating the full, non 
perturbative, field equation for E, still within the sepa- 
rated form (|TT|), that is, assuming the resulting state to 
be still an cigenstate of the angular momentum. The ex- 
istence of the non linear term would then imply \Yi m \ 2 
to be a constant, i.e., independent of both angular vari- 
ables 9 and <j>. Thus, only the case i = can satisfy the 
non linear equation, so that the background a field has 
vanishing angular momentum. 

We now therefore consider again the ansatz ([ll]) with 
no spherical harmonics included and in the region of pa- 
rameter space where a condensate might form, i.e. we 
assume m 2 < 2fr/ 2 . Setting w — yields the actual 
vacuum state as the solution of the field equations 

= ^K 2 + \^ 2 h(h 2 - 1) + 2/^a 2 , (18) 
r z dr dr r z 



1 



q 2 r] 2 Kh 2 + —K(K 2 -1), 



d 2 K 
dr 2 

with boundary conditions for a as 

^(0) = 0, lim o- = 0. 

dr r— too 

Rescalling the field a through 



Y(p) = VA, 



and defining the dimcnsionless parameters as 



Oil 



X a i] 2 



a 2 = 



fm 2 c 



a 3 



AaA^r? 4 



(19) 
(20) 

(21) 
(22) 

(23) 



allows a numerical calculation of the vacuum solution. 
Such a solution, obtained by means of a Successive Over 
Relaxation method |E| , is shown for a special set of pa- 
rameters {oti} on Figure (here and in what follows, the 
background parameter q /X^ has been fixed to the arbi- 
trary value 0.1). 
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FIG. 4. Fields function h(p) [full line], K(p) [dotted line] 
and Y(p) [dashed line] for the set of parameters qi = 0.1, 
«2 = 0.6 and 03 = 0.6. 

These fields represent the vacuum state out of which 
the bound state solutions can be calculated. 



IV. QUANTUM THEORY 

In order to calculate quantum effects related to the 
monopole solution, we turn to the standard solitonic ap- 
proach |2^] . We shall for now on consider the Higgs and 
gauge vector fields to represent a fixed background in 
which the charged scalar field E evolves. This means we 
describes its dynamics through an effective Lagrangian 
density 



£ s = --|^E| 2 -/($ $«-i7 2 )|E| 2 



-^l|E| 2 -— |E| 4 , 



(24) 



2 1 1 4 

and the localized solutions are derivable from the effective 
potential V[E] 

L s = f d 3 x£x = I d 3 x^ - F[E], (25) 



d 3 x 



{i(VE )2 



|E| 2 }. 



(26) 



The classical solution E c derived earlier is then obtained 
by minimizing V[E]: 



5V 
£E 







AE C = [2/($ a $ a - v 2 



m^]E c + A CT E^, 



(27) 
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where, because of the U(l) symmetry in the scalar field 
£, the later could have been chosen real. 

Expanding V[£] around the classical solution and not- 
ing respectively 



and 



V2 



(28) 



(29) 



with ip the quantum perturbation from the classical so- 
lution £ = £ c + ip, one gets 



vm = v[z c \+v in t[fj ) fi\ 



(30) 



d 3 xfj 



/($ Q $ a - V 2 



rug 



where Vi n t[f), p] comprises the interaction terms between 
fj and /I, originating from the self coupling |£| 4 . These 
terms are not explicitely developed here since we are 
only interested in the bound state solutions around the 
monopole, i.e., the stationnary solutions on the basis of 
which the system can be quantized. We are therefore 
looking for the eigenmodes of the second derivatives of 
V with respect to the fields fj and jl, which are then seen 
to satisfy the Schrodinger-like equations 



and 



A + 2/($ a $ a - r? 2 ) 



A + 2/($ a $ a -ry 2 



3E 2 . m 



(31) 



fii = tijN- ( 32 ) 



The eigenmodes fji are of two different kinds: either u>i 
belongs to a discrete set, which is the case if loi < m CT , 
so that fji represents a bound state localized around the 
monopole, or Ui > m a can be parametrized by a continu- 
ous index q that can be identified with a momentum. In 
the later case, the corresponding modes are diffusion, i.e. 
asymptotically free states, with momentum q. Similar 
considerations obviously apply also to the states 

On the basis of this set of eigenmodes, and following 
Goldstone and Jackiw [p4| , one can in principle build 
the quantum theory of the monopoles by constructing 
a Fock space as follows. The original monopole state, 
i.e., the classical solution with no bound state, can be 
boosted to aquire an arbitrary momentum P, thereby 
generating the set {|P)}. The same can be done for 
monopoles with any occupation numbers in the bound 

states {\P,n\ ,rij )}. Finally, the Fock space is com- 
pleted by inclusion of the diffusion states labelled by the 
momenta of the various ingoing and outgoing particles 

{ki^, ' ' ' > <1n > 1i > - ' • >Pm )}• Tllus > an arbitrary quan- 
tum state belongs to the set 



{|iV. 

,(17) 



,(v) 



>) 
l 3 ' 



l 3 



,q N ,<?i , 
eR 3 }, 



(33) 



and this space is assumed disconnected from the ordi- 
nary free particle space of the theory without a monopole. 
Cross sections will then be calculable by computing ma- 
trix elements of the quantum field |S| 4 , which is seen to 
involve essentially space integrals of the eigenfunctions. 
Specific such calculations and their application to primor- 
dial cosmology will be the subject of a further work, and 
for now on we shall concentrate on the actual structure 
of both the self potential and the bound states. 



V. THE SELF POTENTIALS; BOUND STATES 




FIG. 5. Confining potential of the first [F' M ' (r)] and second 
[V {v) (r)} kind in units of the E- field mass squared for 
the set of parameters ai = 0.1, 02 = 0.2 and 03 = 0.1 as 
a function of the rescalled distance to the monopole core p. 
In this case, the minimum of F (M) (r) is located at p = 
while that of V^"(r) is for p = p m 7^ 0. The bound states 
eigenvalues are indicated as the straight lines. The full lines 
correspond to 1/' M ' and its associated eigenvalues, while the 
dotted lines stand for V^ n \ 

The purpose of this section is to exhibit the various 
possible situations, depending on the underlying micro- 
scopic parameters, in which the field £ can evolve. To 
achieve this, we rewrite the eigenmode equations for f) 
and /I in the form 



[-A + V^(r)]fj= {uj 2 ~ 
[-A + VW(r)]/i = (uj 2 - 



ml)fi 



(34) 



which define the self potentials V^(r) and ^'(r) by 
comparison with Eqs. (31) and (p2~|). Some characteris- 



tic shapes and amplitudes of these potentials are shown 
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on Figure pi in rescalled units (the square of the carrier 
mass) , together with the corresponding eigenvalues in the 
same units. 

The figure reveals the existence of two different kinds of 
self potentials, respectively called of the first and second 
kind, depending on the underlying microscopic parame- 
ters. The first kind is characterized by an absolute min- 
imum in the monopole core and corresponds to a weak 
coupling where / is small, while the second kind, hav- 
ing a local maximum at r = and a minimum for some 
nonzero value of the distance to the monopole core, re- 
flects the existence of a strong coupling. As a rule, and 
as could have been expected, the energy eigenstates are 
more bound for confining potentials of the first kind that 
for second kind. Therefore, one can expect the lifetime 
of the corresponding configurations to be strongly depen- 
dent on the parameters. 

The figure |^ illustrates some wavefunctions living in 
self potentials of the first and second kinds. 
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P 

FIG. 6. An illustrative example of some wavefunctions liv- 
ing in a self potential of the second kind. The set of parame- 
ters has been chosen as on figure |^. 

These wavefunctions are all that is needed to fully clar- 
ify the internal structure of the monopole. The expected 
spectrum would have line features as obtained by revers- 
ing the axes on Fig. [|. We therefore turn to these cross- 
sections before discussing the acceleration and propaga- 
tion mechanisms of monopoles. 

VI. MONOPOLE-PROTON INTERACTIONS 

Contrary to the case of an ordinary (structureless) 
monopole, the configurations we have been studying have 
the ability to interact through deep-inelastic scattering 
with air nucleus, a proton say to simplify matters. In- 
deed, a particle trapped inside the monopole core can 



be scattered off the monopole, effectively ionizing it. As 
in the vorton case |pcfl , three possible interactions can 
take place, depending on the proton energy (we keep 
working in the monopole center-of-mass frame): elastic 
interaction, which was discussed in J7J and which will 
presumably yield an unobservable shower accompanied 
with Cerenkov radiation, excitation interaction whereby 
a charge carrier will use the proton energy to move to a 
higher energy level, and finally ionization. The last two 
cases are the most important as one expects them to yield 
an observable signal; they will be compared to existing 
data and used for making predictions in the Pierre Auger 
Observatory Q. 

Once a trapped particle has been moved up to an ex- 
cited state, it will decay into the minimum energy state, 
for instance in radiating a photon^, assuming there does 
exist such a channel. Calling AE — E n — Eq ~ m a the 
energy difference between the two eigenenergies (see, e.g. 
Fig. g), such an observed primary photon would have 
energy, seen on earth 

E obs = 1 AE 1 (35) 

with 7 = s/m M the monopole Lorentz factor. Here, e is 
the total monopole energy, while m M ~ r]/e is its mass. 
It should be remarked that the monopole mass can be 
much larger than the carrier mass, in practice as high 
as a hundred times, so that the monopole itself must be 
accelerated to much higher energies than the observed 
10 20 eV. We shall see however in the following section 
that this is not a very serious constraint. If the energy 
is sufficient to ionize the monopole, the observed energy 
will be also of the order of that given on Eq. (|35|), with 
AE simply equal to Eq, the minimum energy state. In 
order that the reaction actually takes place, it is nec- 
essary that the incident proton, which, as seen from the 
monopole frame, is having the same Lorentz factor 7, has 
enough energy to excite the corresponding state. There- 
fore, AE = 7TO p , with m p the proton mass, so that 

7= ( — ) 1/2 ~ 10 6 - 5 . (36) 
m p 

Such a number would seem to imply states in a monopole 
at energies of the order of a few thousands TeV, and 
therefore far below the cosmological limit m M < 10 9 
GeV. This leaves a large range of possibilities for the 
ratio r}/m a . 

Various cases have to be taken into account if one 
wants to actually understand the data in terms of bound 
states in monopoles. First of all, the determination of 



"The term "photon" shall be used here in the generic sense 
of an unbound decay product of our E— field. In case E is 
effectively coupled to electromagnetism, this "photon" would 
be an actual photon. 
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the energy levels, as discussed above, although specifi- 
cally relevant to our case of a bosonic condensate, is es- 
sentially the same in the case of fermionic bound states. 
However, in the latter case, the Pauli exclusion principle 
would apply and the overall distribution of populated 
states would be different: as all bosons would preferen- 
tially occupy the same state (the lowest energy state), 
fermions on the other hand would fill all the states. The 
expected spectrum would therefore be qualitatively dif- 
ferent. The solid prediction here is the same as that for 
the vorton solution to the UHECR problem j2fj] , namely 
that of the existence of a line spectrum, the exact form 
of which needing a specific model to be determined. If 
such a detection was indeed achieved, the method pre- 
sented above would prove uscfull to actually determine 
the relevant particle physics model through a fit of the 
data. 

Identifying the primary as a monopole is also possible 
through its interactions with air nuclei. In other words, 
one needs to know the monopole atom cross section, or, 
as emphasized above, the monopole proton cross section. 
Such a calculation is feasible in principle by making use 
of the quantum theory developed in Section IV and V, 
and will have to be done in case a ray spectrum is indeed 
observed in the Pierre Auger Observatory H). For the 
time being, a rough evaluation will be enough. 

A monopole such as described here interacts with a 
fermionic field such as a proton through possible exi- 
tations of the Aq modes, the so-called dyonic modes. 
This is a resonant exitation which effectively confines the 
fermion in the core of the monopole for a long time, mak- 
ing the cross-section essentially hadronic (the proton is 
then seen by the monopole as a quasi-bound state) . Once 
trapped, the proton (or a quark therein) has time to in- 
teract with the bound states. As a result [^3|, the inter- 
action cross-section we are seeking is given by: 

% ~ ?~ 2m * 2 > (37) 

numerically of the order of a few hundreds mb. Such 
a large cross-section implies that most of the monopoles 
arriving in the atmosphere would be detectable, provided 
they carry bound states. 



Accelerating a magnetic monopole in the presence of 
a magnetic field is a very easy task: even the galactic 
magnetic field could do it as the kinetic energy E M 
of a monopole in a magnetic field B ~ 10 -6 G with a 
coherence length D ~ 300 pc would be of order 

E M = q,,BL - 6 x 10 19 eV ( ) ( — — ) , 

M Hm \3x io- 6 g; V3oop C y ' 

(38) 

with q M the magnetic charge (inversely proportional to 
the electric coupling constant). Using Eq. (p8[), one gets 
table H for the maximum energy acquired by a magnetic 
monopole (adapted from |25|). 



Magnetic fields, distance and energy 



object 


B 


(G) 


D (pc) 


E M 


(eV) 


Neutron Star 


1Q 9.5 


-10 13 


ltr 12 


10 20 


-10 24 


White Dwarf 


10 4 - 


- 10 s 


ltr 9 


10 18 


- 10 22 


AGN 


10 3 - 


10*5 


10" 3 


10 23 


-10 24 


SNR 


10" 5 - 


- ltr 4 


10 


10 19 


-10 20 


RG lobes 


10" 5 - 


- ltr 4 


10 6 


10 23 


-10 24 



TABLE I. Characteristic values of magnetic fields B and 
associated coherence lengths D for astrophysical objects can- 
didates for accelerating UHECR. The third column gives the 
kinetic energy range that can be obtained for a monopole in 
such a field. 



VII. ACCELERATION AND PROPAGATION 

Once the internal structure of a monopole is known, 
it must be shown that they have the ability to be de- 
tected as UHECR. It is the purpose of this section to 
show that indeed many astrophysical sites have the abil- 
ity to accelerate them. Moreover, there is no GZK cutoff 
for them since they share the vorton features of being 
effectively highly charged, a fact that is compensated 
by their huge mass: the electromagnetic energy losses 
are similar to that of a heavy nucleus having an electric 
charge Z = l/2q = 137/2 j7j, but are reduced by inverse 
powers of the mass pffl . 



7 



Monopoles are however believed not to be bound to 
small objects such as stars (See Ref. @] and references 
therein). Therefore, among the various candidates pre- 
sented on Table |. only Active Galactic Nuclei (AGN) and 
Radio Galaxy (RG) lobes have a chance to be the likely 
accelerators of monopoles. In both cases, the situation 
is such that the regions where the magnetic field is im- 
portant correspond also to regions where other particles 
would be accelerated. These regions are believed to be 
filled essentially with electromagnetic radiation together 
with very low hadronic density (in the sense that protons 
in such a medium would dominantly see the radiation 
and would otherwise have a mean free path, with respect 
to proton-proton interactions, larger than the accelera- 
tion zone). In turn, this radiation is responsible for de- 
grading the accelerated proton energies, but is far below 
the threshold for interacting efficiently with monopoles. 
Therefore, the mean free path of the monopoles is larger 
than the actual size of the acceleration region. As a re- 
sult, one does not expect any kind of cutoff in the injec- 
tion spectrum. 

Finally, one needs to know the expected number of 
such events. Let us assume for the sake of the argument 
that AGN are the most likely accelerator candidates. 
To estimate the flux $ AGN , we set 4> AGN the number 
of monopoles emitted per unit time by a characteristic 
AGN, N AGN the number of AGNs up to a distance which 
we arbitrarily assume corresponds to a redshift z = 1 and 
(D) the mean distance to the AGN under consideration. 
One then obtains 

®ac N = <t>Ac N N AGN (D)- 2 . (39) 

The quantity 4> AGN is 

<P AGN =£to 1 N M , (40) 

with e the ratio between the power used to accelerate 
monopoles and the electromagnetic luminosity, i.e., a 
measure of the efficiency of our mechanism, while to is 
the characteristic time for emitting a monopole. N M is 
the number of monopoles present in the accelerating zone 
of the AGN. Typically N M — an M L 3 AGN , where a is the 
ratio between the hadronic (pb) and monopole (p M ) den- 
sities; n M is the mean number density of monopoles in the 
universe and L AGN the characteristic size of the acceler- 
ation region. Standard monopole formation mechanisms 
give § 

"» = 1,rl '(lo^) 3 < 41 > 

which implies 

N M ~ 10 26 a( ™ M Y ( -^f^Y. (42) 
VlO u Gev/ vicrV^ v ' 

Considering the monopoles to be relativistics so that 
their velocity is essentially that of light, the escape time 



to can be assumed to be the size L AGN , so that Eqs. (E3) 
and (^) combine into 

lo21 Ki^) S (^) 2 («> 

The sphere with 2 = 1 contains a baryonic mass M B = 
5.8 x 10 54 /i _1 g, with h the Hubble constant Ho in units 
of 75 km-s _1 Mpc _1 . The total number of galaxies in 
such a radius is therefore N g = 3 x 10 10 h^ 1 M 1 ' 1 , M w 
being the mass of the galaxy in units of 10 10 solar masses, 
with N AGN approximately a tenth of this value p^j . The 
expected flux $ M is now given by 

$ AGJV ~ 3 x 10 19 e a e M w m 3 n L 2 _ 3 (^T^) 'cm" 2 • s" 1 , 

(44) 

with «6 = a/10 6 , mn = m M /10 11 GeV and L_3 = 
L AGJV /10~ 3 pc. In Eq. (p4|), it should be remarked that 
the efficiency e can exceed unity and in fact presumably 
depends on the energy scale r\. 

VIII. CONCLUSIONS 

Monopoles represent the most generic prediction of 
GUT models and as such are always presented as a possi- 
ble cosmological nuisance. This is because they are usu- 
ally postulated to form at the GUT phase transition, in 
which case their remnant density in the universe would be 
proportional to the fourth power of this energy scale [|| , 
a density considerably larger than the critical density to- 
day. Going back to the original idea of Kephart and 
Weiler @, we consider instead their usefulness in ex- 
plaining the mystery of the UHECRs by assuming the 
simplest of all solution to the monopole problem, namely 
that they are produced at a phase transition taking place 
at a temperature scale no higher than 10 9 GeV, implying 
a monopole mass m M < 10 11 GeV. In this case, one can 
safely consider that magnetic monopoles do exist in rea- 
sonnable number in the universe and their study becomes 
less academic. 

Monopoles as such can hardly interact with air pro- 
tons to yield air showers as they must be topologically 
stable: their expected signal would be an extremely dif- 
ficult to observe Cerenkov shower [Q. However, they are 
obviously easily accelerated to energies much higher even 
than the world record observed until now of E w R ~ 
3 x 10 20 eV J2(|. One is therefore tempted to consider 
them as candidates for explaining those data. The point 
we want to make here is that a possibility is left opened 
when one considers the nature of monopole. 

A monopole is a solitonic configuration of a winding 
localised Higgs field. In most reasonnable theories, such 
a Higgs field not only permit the symmetry breaking, and 
thus the appearance of topological defects, but serves 
also as a means to provide masses to various particles 
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to which they couple. This in turn implies that these 
particles might get trapped in the monopole core, form- 
ing bound states that might interact with the air pro- 
tons in such a way as to be effectively expelled from the 
monopole. Moreover, they behave as a neutron in a nu- 
cleus: they are stable in the form of bound states, but un- 
stable otherwise. The observations can then be explained 
in supposing that monopoles regularly hit the earth at- 
mosphere with tremendous energies, releasing only part 
of it in such "ionising" processes. The unstable parti- 
cles thus obtained are relativistic, although not highly 
with a Lorentz fact not exceeding 10 7 and can initiate 
air showers. 

We have developped a model in which, to make things 
simple, the bound states are formed by means of the con- 
densation of a scalar field. Such a field can be charged, 
a point that would modify our analysis by corrections of 
order e ~ 1/137, but it should be reminded that the 
leading effects we studied come from the topological de- 
fect itself and are thus of order e -2 ~ 137; the correction 
we neglected is therefore some four orders of magnitude 
smaller. Given this approximation, we have examined in 
detail the internal structure of the charged monopole and 
derived the expected energy levels that might hopefully 
be measured in a precise UHECR experiment, as the lat- 
ter should, in our model, yield a line spectrum. If such an 
observation was done, the formalism we have presented 
would allow a straightforward computation of the free 
parameters and a reconstruction of the underlying the- 
ory. Whether or not the Pierre Auger Observatory ]2lJ| 
project will fulfill this task is yet an open question and 
depends essentially on the spacing of the energy levels; 
for some cases, one expects that it will. 

Our main result is that such a model satisfies all the 
present observational constraints: their expected flux 
is of the correct order of magnitude and they inter- 
act strongly with the atmosphere. This last feature 
comes from the pre sence of bound states and the Callan- 
Rubakov effect |13|] , A model based on such topological 
defects would have the advantages of the non-acceleration 
scenarios (bottom-up models) since they would easily 
propagate, together with the advantages of acceleration 
mechanisms, as they need point sources in order to be 
observed. Note that this is compatible with the most re- 
cent data (observed doublets and triplets of events within 
2.5°) implying localized acceleration. These very 
data render the previous mechanism for UHECR using 
monopoles tt very unprobable. 
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